Accelerating waves in polar coronal holes as seen by EIS and SUMER 

G. R. Gupta^'^, D. Banerjee', L. Teriaca^ S. Imada"^, S. Solanki-''^ 
girjesh@iiap.res. in 

^H . 

^ ■ ABSTRACT 

We present EIS/Hinode & SUMER/SoHO observations of propagating disturbances 

detected in coronal lines in inter-plume and plume regions of a polar coronal hole. The 

^ . observation was carried out on 13''' November 2007 as JOP196/HOP045 programme. The 

SUMER spectroscopic observation gives the information about the fluctuation in radiance 

and on both resolved (Doppler shift) and unresolved (Doppler width) line-of-sight velocities 

^ I whereas EIS 40" wide slot images detect fluctuations only in radiance but maximizes the 

(-H ' probability of overlapping field of view between the two instruments. From distance-time 

Qh' radiance maps, we detect the presence of propagating waves in a polar inter-plume region with 

P ■ a period of 15 min to 20 min and a propagation speed increasing from (130 ± 14) km s~^ just 

above the limb, to (330 ± 140) km s~^ around 160" above the limb. These waves can be traced 

C^ . to originate from a bright region of the on-disk part of the coronal hole where the propagation 

speed is in the range of (25 ± 1.3) km s"^ to (38 ± 4.5) km s"^, with the same periodicity. 

K^ , These on-disk bright regions can be visualized as the base of the coronal funnels. The adjacent 

(^ I plume region also shows the presence of propagating disturbance with the same range of 

^ ' periodicity but with propagation speeds in the range of (135 ± 18) km s~^ to (165 ± 43) km s~^ 

ff^ ■ only. A comparison between the distance-time radiance map of both regions, indicate that the 

l/^ ■ waves within the plumes are not observable (may be getting dissipated) far off-limb whereas 



P^ 



o 
o 



this is not the case in the inter-plume region. A correlation analysis was also performed to find 
out the time delay between the oscillations at several heights in the off-limb region, finding 
results consistent with those from the analysis of the distance-time maps. To our knowledge, 
this result provides first spectroscopic evidence of acceleration of propagating disturbances in 
\-i ' the polar region close to the Sun (within 1.2 R/Rq), which provides clues to the understanding 

of the origin of these waves. We suggest that the waves are likely either Alfvenic or fast mag- 
netoacoustic in the inter-plume and slow magnetoacoustic in plume regions. This may lead to 
the conclusion that inter-plumes are preferred channel for the acceleration of the fast solar wind. 
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1. Introduction 



C oronal holes are regions o f cool and low density plasma that, as such, are 'dark' at coronal temper- 
atures (IMunro & Withbroelll972h . During solar minima, coronal holes are generally confined to the Sun's 
polar regions, while at solar maxima they can also be found at lower latitudes, usually associated with rem- 
nant active regions, as so-called 'equatorial' coronal holes. The pred ominantly unipolar magnetic field from 
coronal hole regions is thought to give rise to the fast solar wind (e.g. jKrieger et alll973h . During solar min- 



imum, Ulysses observations clearly show that the solar wind exhibits two modes of outflow: the fast wind, 
associated with polar coronal holes, with outflow speeds o f » 800 km s~' and the slow wind with outflow 



speeds of ~ 400 km s~' associated with equatorial regions (|Woch et al.lll997l : lMcComas et al.ll2000h . How 
ever, during sol ar maxima, low la titude coronal holes also show faster than average solar wind speed upto 
w 600 km s"^ (|Zhang et alJl2003h . Extreme-ultraviolet im ages of polar coronal holes reveal the presenc e 
of diffuse, spike-like or sheet-like structures called plumes (|Bohlin et alJI 19751 : lAhmad & Withbrodll977n . 
which subtend a n angle of roughly 2° relative to Sun center at low altitude and expands super-radially with 
the coronal hole (IDeForest et al.ll 19970 . Regions between these structures are termed as inter-plumes. From 
VUV spectroscopy, plumes are known to be denser and cooler than the surrounding inter -plume regions 



(e.g...Wilh elm 2006), while spe ctral lines are o bserved to be broader in inter-plumes (i.e., iBanerjee et al. 



2000bl : iGio rdano et al. 2000; Te riaca et al.ll2003h . However differences in mass, momentum and energy flux 
in plumes and in inter-plumes are still not known precisely. 

There are several theoretical models which describe the role of MHD waves in the acceleration of the 



fast solar wind in coronal holes (see review by lQfmanll2005l : ICranmeii l2009l and references therein) and 
inter-plumes are often believed to be the primary site for thi s acceleration. It is further conjectured that 



these waves originate from the on-disk brig ht network regions dWang et al. 1997 1 Pats ourakos & ViafcOOO ; 



Giordano et al. 



2000 



Banerjee et al.ll200lh . A number of studies d Ofman et al.ll 19971. i2000i : iBanerjee et al. 



2001, : Popescu et al.l 12005.) have reported detection of oscillations in the off-limb regions of polar coronal 
holes. All of t hese studies point to the presence of compressional waves, thought to be slow magneto- 
acoustic waves (IDeForest & Gurmanlll998l : lo'Sheaet allbooel 120071 : Banerjee et all 2009bl). On th e other 
hand, e vidence for Alfven wave s propagating into the cor ona had been reported by IBanerjee et all (|l998 . 
2009ah : lDolla & Solomonl(|2008h : lLandi & Cranmen (|2009r) by studying the line width variations with height 
in polar coronal holes. Rece nt reports of detection s of low-frequency (< 5 mHz), propagating transverse 
motions in the solar cor ona (ITomczyk et all l2007h (from coronagraphic observation) and chromosphere 



(DePontieu et al. 



(IDe Pontieu et al. 



2007bl) and their relationship with chromospheric spi cules observed at the solar limb 
2007 ah with t he Solar Optical T elescope aboard Hinode (IKosugi et al.l2007h have widened 
aul 



interest in the subject. Recently iJess et al.l (120091) have reported detection of torsional Alfvenic motions as- 
sociated with a large on-disk bright-point group. These wav es are believed to be a promising candid ate for 
the heating of the corona and acceleration of the solar wind (|Belcherlll97ll : ISuzuki & Inutsukall2005h . 



Furthermore, it has been suggested that the fast solar w ind streams ori ginate from coronal hole funnels 
and are launched by reconnection at network boundaries, (ITu et all 120050 . Measurements of the outflow 
speed in the extend ed corona have been obtained with the Ult r aviolet Coronagraph Spectro meter (UVCS ) 
aboard SoHO (e.g., Untonucci etallboool : Ixeriacaet allboOBi : Untonucci etal] 12004. : .Telloni et allboOTh 



Some of these studies conc l uded that plurnes have lower outflow speeds than inter-plume regions (INoci et al 



1997l : lGiordano et alboool i lwilhehn et allboool i lpatsourakos & ViaJ20od : lTeriaca et alboOsl i lRaouafi et aL 



2007h and, hence, may not contribute significantly to the fast solar wind, whereas some other theoretical and 
observational studies find hig her outflow speeds in plumes than in inter-plume regions for at least some alti- 



oDservational stuaies nna mg ner outilow speeds m ptumes man m mter-ptume regions tor at least some alti- 
tudes above the photosphere dCasalbuoni et al.ll999l : lGabriel et al.ll2003ll2005r) . These contradictory reports 



led to the debate on whether plumes or inter-plumes are the preferred source regions for the acceleration of 
the fast solar wind. This topic is highly debated and still open for further confirmation. 



Recently, iBanerjee et al.l (|2009bh reported the detection of propagating slow magnetoacoustic waves 
with periods between 10 min and 30 min and speed ~ 75 km s"' to 125 km s"^ above the limb of a polar 
coronal hole. In their study, the propagating disturbances which are due to radiance perturbations are seen 
from the limb region up to « 100" above the limb. There is no discernible acceleration or deceleration of 
any individual feature as it propagates. In that study, the oscillations were detected in the two spectral lines 
of Ne vm 770 A and Fe xii 195 A observed with the Solar Ultraviolet Measurements Of Emitted Radiation 



(SUMER. IWilhelm et al.ll 19951) aboard the Sola r and Heliospheric Observatory (SoHO) and with the EUV 
imaging spectrometer (ElS. ICulhane et al.ll2007r) aboard Hinode, respectively. 



In this paper, we combine again the capabilities of SUMER and EIS to observe the on-disk, limb and 
far off-limb region of the coronal hole, to search for the origin of waves close to the Sun and study their 
propagating nature. The plan of the paper is as follows: in section [2l the observations acquired for this study 
and the data reduction techniques are outlined. In section [3] results of the present study are presented with 
the distance-time radiance map analysis, power series analysis and time-delay analysis. A discussion of the 
observational results and a comparison with similar results are taken up in section |4] and finally conclusions 
are drawn in section [5] 



2. Observations 



2.1. Data 

The data analyzed here were obtained on 13*'' November 2007 during a Hinode/SUMER joint observing 
campaign as part of the Hinode Observing Programme (HOP) 45/Joint Observing program (JOP) 196. The 
data consist of time series taken by SUMER and EIS in the north polar coronal hole. For SUMER, the 
1" X 120" slit was centered on the limb and spectral profiles of the Ne vni 770 A, O iv 790 A and S v 786 A 
were acquired from 19:13 to 22:15 UTC with an average cadence of 18.12 s in sit and stare mode. The 
exposure time was 18 s and a total of 600 time-frames were obtained during the observation. For EIS, the 
40" wide slot was used to obtain 40" x 5 12" images in several spectral lines in the wavelength ranges of 
170 to 210 A and 250 to 290 A with spatial resolution of 1" pixel"^ over the time interval from 18:20 to 



23:50 UTC. The exposure time was 45 s with an effective cadence of w 47 s. A total of 420 time-frames 
were obtained during the observation. Before the start of the temporal series, raster images were obtained 
with SUMER and EIS in order to co-align and to provide context. During the observation, the EIS slot 
covered the quiet Sun south of the coronal hole as well as the on-disk and off-limb parts of the hole (see 
Figure [T]). Table [J lists the emission lines included in this study from EIS and SUMER, their formation 
temperature, and the location where the radiance maximum is observed. Top left panel of Figure [T] shows 
the location of the different slits on an FIT image taken in the Fe xii 195 A passband. The rectangular box 
marks the location of the EIS slot while the dashed line gives the location of the SUMER slit. The radiance 
variation along the solar-X at solar-Yw 1000" is over-plotted as a white line in arbitrary units and allows us 
to identify the locations of plume and inter-plume regions within our field of view revealing that the SUMER 
slit is pointed within an inter-plume region while the EIS slot covers both plume and inter-plume regions. 
Bottom left shows the context raster obtained by EIS in Fe xu line whereas right panels show the context 
rasters in O iv (top) and Ne vm (bottom ) spectral lines as o btained by SUMER. Figure |2] corresponds to 



the images obtained by EUVI/STEREO (IHoward et al.ll2008h for the same region. The angular separation 



between the two spacecraft, about 40 °, allows an estimate of the orientation of the plume. 

Table 1 : Emission lines observed with EIS and SUMER and position of the respective limb brightening 



Ion 


Wavelength (A) 


log Ti^ax (K) 


Limb brightening 


Hen 


256.32 


4.9 


985" 


Sv 


786.47 


5.2 


988" 


Oiv 


790.19 


5.2 


989" 


Mg VI 


270.39 


5.6 


990" 


Fe vni 


185.21 


5.6 


990" 


Mg vu 


278.40 


5.8 


991" 


Sivn 


275.35 


5.8 


993" 


Ne vm 


770.42 


5.8 


993" 


Fex 


190.04 


6.0 


993" 


Fexi 


188.23 


6.1 


- 


Fe xu 


195.12 


6.1 


992" 



2.2. Data reduction and alignment 

All data have been reduced and calibrated with the standard procedures given in the SolarSoft (SSWJj 
library. SUMER data were first decompressed, corrected for response inhomogeneities (fiatfield), dead- 
time, local-gain and for geometrical dis tortion (de-stretch), using the most recent standard routines (see 



WiUielm et al.lll997l : iTeriaca et al.lll999l) . After these steps, data still showed a residual pattern from the 



micro-channel plate structure that was removed using a correction matrix obtained by first averaging all 



Ihttp : //sohowww . nascom . nasa . gov/solar sof t/ | 



spectral images and then applying a low pass filter to the average. Single Gaussian fitting was used to 
retrieve the line amplitude, position and line widths of the SUMER spectral lines. Before fitting, a running 
average over three pixels along the slit and over three consecutive spectra was applied to improve the signal 
to noise ratio of the SUMER data. Line positions from the fitting are then converted into Doppler shifts by 
taking as reference the average over the disk part of the image. The observed line widths are corrected for the 
instrumental profile by applying a de-convolution function taking into account the order of diffraction and 
the slit width used during the observations (SolarSoft routine con-width-funct-A). The raw EIS data were 
processed by the standard SolarSoft programme eis-prep which helps in removing detector bias and dark 
current, hot pixels and cosmic rays and returns absolutely calibrated data. The movement of the slot image 
on the detector due to thermal variations along the orbit was conected. The displacement in the dispersion 
(solar X) direction was obtained by measuring the position of the edge of the Fe xu 195 A slot ima ge ove r 
time. The displacement in the Y direction is taken equal to 2.5 times that in the X direction dlmadalbOlCl) . 
The validity of the latter assumption was verified by checking the limb Y position vs. time. Finally, EIS 
data were corrected for the spacecraft jitter by using housekeeping data. Figure [3] shows the variation of 
the radiance along a vertical strip (about 100 pixels long, centered over the limb) at the overlap position of 
the EIS slot with the SUMER sht as a function of time (x-t slice). The x-t slice are for data without any 
correction (top panel), with the orbital effects corrected (middle), and with correction for both orbital and 
jitter effects (bottom). Data are also affected by passage of the spacecraft through South Atlantic Anomaly 
(SAA). The regions of affected data appear around 20:35 and 22:15 UTC and, less clearly visible, around 
18:55 UTC. These affected portions are replaced by linear interpolation. It should be finally noticed that 
for wavelet analysis, EIS data have been binned over 5 and 9 pixels in the X and Y-directions, respectively. 
This will further smooths out any residual jitter or orbital variation. The data from the short wavelength 
detector were shifted in the Y-direction to compensate the wavelength dependent offset between the short 
and long wavelength detectors (Kamio & Hara, private communication). To align the different instruments, 
the SUMER context raster has been chosen as the reference. Hence, the EIS context raster has been cross- 
correlated with the SUMER raster and an offset of 9" in the E-W direction (solar X), and -24" in the N-S 
direction (solar Y) was found and corrected for. After the alignment, the pointing of the different instruments 
are plotted in Figure [T] We have plotted the variation of radiance along the solar-Y at the overlap region of 
the SUMER slit and EIS slot at solar-X « -72", see Figured This allows us to identify the location of the 
limb brightening in various spectral lines and their radiance fall-off in the off-limb region. We set the limb 
position at solar- Yw 985", as identified from the limb brightening of He ii in Figure |4l The radiance peaks 
of the He u line for the disk part of the coronal hole are used to identify bright locations (presumably the 
footpoint of the coronal funnels). 



3. Results 

In this section we will present detailed analysis of the sit-and-stare observation as recorded by EIS and 
SUMER. From the context image. Figure [TJ it is clear that the SUMER slit is pointing to an inter-plume 
region and the EIS slot covers both plume and inter-plume region, overlapping with the SUMER slit. Thus, 
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Fig. 1. — Top left: The location of the different slits are over-plotted on the EIT image taken on 13*'' 
November 2007 in the passband dominated by the Fe xii 195 A line. The rectangular box marks the location 
of the EIS slot while the dashed line gives the location of the SUMER slit. The EIT radiance variation along 
the solar-X is over-plotted as a white line in arbitrary units at fixed solar-Yw 1000". This variation along the 
solar-X allows us to identify the location of plume and inter-plume regions. Bottom left: The EIS context 
raster taken in the same Fe xii line shows the location of both the EIS slot (rectangular box) and the SUMER 
slit (dashed line) during the sit-and-stare sequence. Right: The SUMER context rasters taken in the O iv 
(top) and Ne vui (bottom) spectral lines. The continuous line gives the radiance variation along the solar-X 
in an arbitrary unit. In all panels, the arrow indicates the location of the bright region from where waves are 
presumably originating. 
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Fig. 2.— EUVI/STEREO and EIT/SoHO images in the 171 Abandpass taken around 19:00 UT on U"' 
November 2007. The solid line shows the plume axis as obtained by fitting the positions of the radiance 
maxima of horizontal cuts on the EIT image. Separation between the two STEREO spacecraft was 40° (20° 
with Earth/SoHO). The diff'erent angular scale (km/") was taken into consideration. Numbers from 1 to 4 
on the EIT image identify reference elements that can be recognized on the EUVI-B image (suffix b) and/or 
on the EUVI-A image (suffix A). 

in the inter-plume region, EIS and SUMER will both provide information whereas in the plume region only 
EIS observations are available. Analysis of the radiance maps (x-t slices), wavelet analysis and correlation 
analysis will be performed at the two locations as described in the following three subsections. 



3.1. Radiance x-t slices 



Maps of radiance along the slit vs. time (distance-time map or x-t slices) were built using the SUMER 
Ne vm integrated line radiance and the EIS Fe xu radiance averaged over 5" in the X direction at the po- 
sition overlapping with the SUMER slit (which covers the inter-plume region). The resulting maps were 
then smoothed over « 3 min and the background trend of w 20 min has been subtracted from each solar-Y 
pixel along time. In general, similar procedures are applied while doing the Fourier or wavelet analysis of 
a time series. In the x-t slices, the presence of alternate bright and dark regions indicate the presence of os- 
cillations. Moreover, diagonal or slanted radiance enhancements are signature of propagating disturbances. 
Thus, from such maps it b ecomes possible to estimate periods and projected propagation speeds (see e.g. 
DeForest & GurmaJl998l) . 



Both the SUMER slit and the EIS slot are centered on the solar limb and, hence, cover the region 
on-disk as well as off-limb. As the observed region is near to the pole, the effect of solar rotation is very 
small and amount to less than about 3" per hour at 100" below the limb. We first concentrate our attention 
to the inter-plume location, around solar-Xw -72", probed by both instruments. Figure [5] shows the x-t map 
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Fig. 3. — Variation of radiance over time along the slice of corona recorded by both SUMER and EIS. The 
top panel shows the distance-time map of the Fe xn Une as recorded by EIS. Here the limb clearly shows 
the effect due to orbital variation and instrumental jitter. The middle panel shows the distance-time map 
after orbit correction and the bottom panel shows the map after orbit and jitter correction applied, showing 
the limb position to be much more stable in solar-Y as compared to the uncorrected data. In the panels, 
two visible vertical stripes around 20:35 and 22:15 UTC are damaged data, probably due to SAA transits 
(another, less visible stripe is present around 18:55 UTC). 
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Fig. 4. — Time averaged radiance variations along the SUMER slit and EIS slot (solar- Y) at solar-X- 
for different Lines as labeled. 
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of the radiance of the Ne viii spectral line, where the presence of slanted bright and dark region is clearly 
visible. A disturbance appears from the on-disk bright region around solar-Y « 967" (see Figure [T] and |4ll, 
and propagates towards the limb. The signature of oscillations ai^e very strong in this bright region. No 
signature of propagation is visible below solar-Yw 967" (the SUMER slit covers down to solar-Yw 930"). 
Hence the assumption that this bright region is the source of these propagating disturbances is justified 
(Figure [5]l. The speed of propagation measured from the slope of the enhanced slanted radiance stripes is 
(25 ± 1.3) km s'^. This average speed is measured up to solar-Yw 992" which is very close to the region 
of limb brightening (see also Figure |4l). As the propagation reaches the limb, its speed increases and the 
enhanced features become more vertical. Up to solar-Yw 1010", the measured speed is (38 ± 4.5) km s"'. 
This change in speed is a clear signature of acceleration of the propagating disturbance. Furthermore, when 
this propagation reaches beyond the limb brightening height i.e. in the corona, its speed further increases 
to (130 ± 51) km s~' and the propagation is seen up to solar-Yw 1020". Beyond this height the signature 
of propagation becomes very poor, most likely due to the low signal. The periodicity of the fluctuations is 
w 14 min to 20 min. In Figure [51 the over-plotted white lines follow the slope of the enhancements and 
are plotted with a periodicity of w 14 min. It can be seen that in some places the over-plotted white line 
do not coincide with the enhanced lanes but it is nevertheless parallel to it. This suggests that even if the 
periodicity changes within a certain range, the propagation speeds are fairly uniform. There is no clear 
evidence of propagating disturbances (in terms of velocity fluctuations) in the line-of-sight (LOS) velocity 
(obtained from Doppler shift) x-t map, although periodicities are revealed by the wavelet analysis (see next 
section). 

Figure [6] shows the radiance x-t map of the Fe xn 195 A line at the same location (solar-X ~ -72", 
inter-plume region). The image was processed as described above. The analysis of the x-t map over the 
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Fig. 5. — Enhanced distance-time (x-t) map of radiance (along solar-Y) variation at solar Xw -72" as 
recorded by SUMER in the Ne vm spectral line on 13''' November 2007. Here the slit covers the on-disk, 
limb and off-limb region of the polar coronal hole and it is positioned in the inter-plume region. The slanted 
lines correspond to the disturbances propagating outward with increasing speed. The dashed horizontal line 
indicates the position of the limb brightening in Ne vm. In the on-disk region the disturbance propagates 
with a speed of (25 ± 1.3) km s~\ increasing to (38 ± 4.5) km s~^ close to the limb and to about 
(130 ±51) km s~' in the off-limb region. The periodicity is in the range of w 14 min to 20 min as also 
obtained from the wavelet analysis, see Figures [8] & |9] 

coronal hole, reveals that there is no clear signature of propagating disturbances in the on-disk coronal hole 
region. These alternate bright and dark regions are clearly visible only around the limb and far off-limb, 
hence only these regions are plotted in Figure [6l In this map, the propagating disturbances are visible from 
solar- Yw 1000" up to the upper end of the EIS slot at solar- Y« 1 140". In this inter-plume region, the dis- 
turbance propagates with an average speed of (130 + 14) km s~' from solar-Y« 1000" to solar-Yw 1085" 
increasing to an average speed of (330 ± 140) km s~^ up to solar-Y« 1 135", clearly showing signatures of 
acceleration. The periodicity of the fluctuations is in the range of « 15 min to 18 min. Here the over-plotted 
white lines gives the slope of the enhancements and are plotted with a periodicity of « 17 min. It can be 
again seen here that in some places the over-plotted white line do not coincide with the enhanced lanes but 
is nevertheless parallel to it as was discussed for the results from the Ne vm Une. 



Summarizing, observations in the Ne vm spectral line reveal a propagating disturbance originating 
in a bright region in the on-disk coronal hole that starts propagating towards the limb region with speed 
(25 ± 1.3) km s~^ Near the limb the speed increases to (38 ± 4.5) km s"^ and reaches (130 ±51) km s~^ in 
the off-limb region. A similar speed is measured at the same height in the Fe xn spectral line by EIS. Hence, 
both instruments see approximatively the same speed in the same region, which is different from the result 



11 



1140 
USD 




EIS - 


FeXIE 














HH 


1100 










pPBUi 


J^ loeo 




^hSE!''' 




^fc.-?.^^B 


t* 


-1 1000 —^^K^^I^E 


£^:l 






pjBl 


1040 '~^^^nKJ^^^'' 


5k •< 






«sS 




- ^^-i^^m^a 


^^E-j 






k. ^^^1 


lOEO 
lODO 


@ei 


Pi 









lfi:DO 1B;S0 

Time (UT) 




S3;40 S^rQO 

Itaie (UT) 



Fig. 6. — Enhanced x-t map of radiance variation along solar-Y at solar-Xw -72" as recorded by EIS in 
Fe XII on 13^'' November 2007. The height range shown here covers the near off-limb and far off-limb 
region of the polar coronal hole and corresponds to the inter-plume region. The slanted hnes correspond 
to the disturbances propagating outward with increasing speed. In the near off-limb region the disturbance 
propagates with speed of (130 ± 14) km s"^, and accelerates to (330 ± 140) km s"^ in the far off-limb 
region. The periodicity is in the range of w 15 min to 18 min as obtained from the wavelet analysis, see 
Figure nn 
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reported by iBanerjee et alj (l2009bl) for a likely plume region, where they find different propagation speeds 



in different lines. Further off-limb, the speed of the propagating disturbance reaches (330 ± 140) km s~' as 
seen by EIS in the Fe xu spectral line. Overall, acceleration of propagating disturbances are observed from 
on-disk to far off-limb in an inter-plume region simultaneously by two different instruments on different 
satellites. 
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Fig. 7. — Enhanced x-t map of radiance variation along solar-Y at solar-Xw -39" as recorded by EIS in 
Fe XII on 13'^^ November 2007. The height range shown here covers the near off-limb and far off-limb region 
of the polar coronal hole and falls in the plume region. The slanted lines corresponds to the disturbances 
propagating outward with nearly constant speed. In the near off-limb region the disturbance propagates with 
speed of (135 ± 18) km s~', and accelerates to (165 ± 43)kms"^ in the far off-limb region. The periodicity 
is in the range of w 15 min to 20 min as obtained from the wavelet analysis, see Figure [ 



The EIS slot also covers part of a plume and the position around solar-Xw -39" is selected to represent 
this region. The processed x-t radiance map is plotted in Figure |7] Also in this case, there is no clear 
signature of propagating disturbances in the on-disk coronal hole region. These alternate bright and dark 
regions are clearly visible only at the limb and off-limb, hence again only these regions are plotted in 
Figure |7] In this map, the propagating disturbances are visible from solar-Yw 1000" up to solar-Yw 1120". 
In this plume region the disturbance propagates with a speed of (135 ±18) km s~' from solar-Yw 1000" to 



-13- 

solar-Yw 1075" and with (165 ± 43) km s~' up to solar-Yw 1120". Beyond this height the map becomes 
diffuse and there is no clear signature of propagation. Although this may be interpretated in terms of wave 
dissipation, this lack of signature at greater heights may be simply due to merging with the background 
signal. 



3.2. Analysis of oscillations 

In this subsection, in order to study the detailed properties of the propagating disturbances as seen 
in the enhanced radiance x-t maps (Figures [51 [6] and IT]), we make use of wavelet analysis and focus on 
individual locations in the on-disk and off-limb corona. The full SUMER time series is used to detect 
oscillations in the radiance, Doppler shift as well as in the Doppler width of the Ne vin spectral hne at several 
locations. The bright location is identified on-disk using the maximum radiance seen in He ii 256 A by EIS. 
At several off-limb locations there is sufficient signal to noise to detect oscillations with a high confidence 
level. Figures [8l |9l [10] & \TT\ show examples of oscillations measured in the polar region at fixed solar- 
X~ -72" (which corresponds to the inter-plume region) and at several solar-Y locations: « 967" (on-disk), 
w 1020" (off-limb, but close to the limb) and « 1120" (far off-limb), as mentioned in the figure caption. On 
the other hand. Figure [12] shows oscillation measured at solar-X~ -39" (which corresponds to the plume 
region) and solar-Yw 1030" (off-limb). In these figures the top panel shows the variation of the radiance 
(hereafter the term radiance will be used for trend subtracted integrated line radiance) with time. Details 
on the wavelet analysis, w hich provides information on the temporal variation of a signal, are described in 



Torrence & Compol (119981 ). For the convolution with the time series in the wavelet transform, the Morlet 
function is choosen. The oscillations shown in the upper panel had their background trend removed by 
subtracting from the original time series a 100-point {~ 30 min) and 35-point (w 30 min) running average 
for SUMER and EIS data, respectively. In the wavelet spectrum, the cross-hatched regions are locations 
where estimates of oscillation period become unreliable which is called as cone-of-influence (COI). As 
a result of the COI, the maximum measurable period is shown by a horizontal dashed line in the global 
wavelet plots, which are obtained by taking the mean over the wavelet time domain. This global wavelet 
is very similar to the Fourier transform as both are giving the distribution of power with respect to period 
or frequency. Whenever the Fourier spectrum is smoothed, it approaches the global wavelet spectrum. The 
period at the location of the maximum in the global wavelet spectrum is printed above the global wavelet 
spectrum. 

From the x-t map analysis (Figure [5]) it was seen that outward propagating disturbances at the inter- 
plume location (solar-Xw -72") originate from a bright on-disk region at solar-Yw 967". As this region 
has been covered by the SUMER slit, we have information about the radiance, Doppler shift as well as the 
Doppler width of the Ne vm spectral line. Time series have been obtained at this bright location by taking 
a 5" average over solar-Y and then wavelet power spectra have been plotted for both radiance and LOS 
velocity in Figure [8] There is a clear presence of « 15 min to 20 min periodicity in both radiance and 
velocity. Going to the off-limb inter-plume at solar- Yw 1020", time series in both radiance and Doppler 
velocity were obtained by averaging over 9" in the Y direction (to increase the signal to noise ratio). Also 
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Fig. 8. — The wavelet result for the on-disk location at solar-Yw 967" and solar-Xw -72" in Ne vm radiance 
(left side) and velocity (right side). In each set, the top panels show the relative (background trend removed) 
radiance/velocity smoothed over 3 min. Bottom left panels show the color inverted wavelet power spectrum 
with 99 % confidence level contours while bottom right panels show the global (averaged over time) wavelet 
power spectrum with 99 % global confidence level drawn. The period PI at the location of the maximum in 
the global wavelet spectrum is printed above the global wavelet spectrum. 
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Fig. 9. — The wavelet analysis results corresponding to solar- Yw 1020" in the Ne vin radiance (left side) and 
in velocity (right side) at solar-Xw -72" (inter-plume region). See the caption of Figure |8]for a description 
of the different panels. 
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Fig. 10. — The wavelet analysis results for the oscillations in Doppler width of the Ne viii line at solar- 
Yw 967" (left side) and at solar-Yw 1020" (right side) obtained at solar-Xw -72" (inter-plume region). See 
the caption of Figure [8] for a description of the different panels. 

from the wavelet power spectra (Figure |9l), a clear presence of « 15 min to 20 min periodicity in both 
radiance and doppler shift was found. Wavelet power spectra of the Doppler width time series were also 
obtained at the two locations and are shown in Figure [TO] The analysis reveals a periodicity similar to that 
observed in radiance and in Doppler shift. 

EIS Fe xu time series have been produced by averaging over 5" and 9" in the X and Y-direction, 
respectively. Then a wavelet analysis was performed at two inter-plume locations; one near the limb (solar- 
Y~ 1020") and the other further off-limb (solar-Yw 1120"). The wavelet power spectra are shown in 
Figure [TT] Both heights show periodicity between 15 min and 20 min, consistent with the results from the 
x-t map (see Figure O. Furthermore, these periods are consistent with the periods obtained from Ne viii. 
In summary, again it can be concluded that the propagating disturbance, which originates from the on-disk 
bright region, propagates off-limb and in the far off-limb inter-plume region. 

In order to check the periods of propagating disturbances in the plume region, the wavelet analysis has 
been carried out also at location solar-Yw 1030" and solar-Xw -39" by averaging over 9" in solar-Y and 
5" in solar-X, (see Figure [T2)) . Also in this case, the period of propagation, wl2 min to 20 min, is consistent 
with the x-t map in Figure |7] However, as seen from the x-t map, these disturbances are not visible at greater 
heights above the limb. 



3.3. Correlation Analysis 

In the earlier subsections we presented results from the Ne vin and Fe xii lines only. Propagation prop- 
erties of waves can also be studied by conelation analysis. First, we focus our attention on the on-disk bright 
region, where we expect the waves seen at the inter-plume location to originate. At this location (Yw 967"), 
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Fig. 11. — The wavelet analysis results corresponding to the Fe xii radiance at solar-Yw 1020" (left side) 
and at solar-Yw 1120" (right side) at solar-Xw -72" (inter-plume region). See the caption of Figure [8] for a 
description of the different panels. 
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Fig. 12. — The wavelet analysis results corresponding to the Fe xii radiance at solar-Yw 1030" and solar- 
Xw -39" (plume region). See the caption of Figure [8] for a description of the different panels. In this case, 
confidence contours are drawn at 95 % level. 
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we find clear presence of oscillations in different lines as recorded by SUMER and EIS (as tabulated in 
Table. [T]l. For the on-disk study, we will concentrate on the correlation between different lines as recorded 
by SUMER and EIS while for the off-limb study we will calculate correlation coefficients between different 
heights as recorded by the same Fe xii line. The correlations between the time series from two different lines 
have been obtained using the IDL routine C-CORRELATE at different time delays between the two series. 
Correlation coefficients have been calculated for six line pairs and are plotted in the left panel of Figure [T3l 
The time resolution is about 18 s for SUMER and 47 s for EIS (governed by the respective cadences). The 
time delay at peak of correlation can be considered as the time delay between the oscillations in the two lines 
forming the pair. It can be seen that the level of correlation and the time delay are in inverse proportion, 
which means that the correlation is higher for the lines having smaller temperature separation and, because 
in coronal hole the temperature changes mainly as a function of height, smaller time delay. The correlation 
between He ii and Fe xii is smaller and shows the largest time delay whereas the correlation between He u 
and Fe x is comparatively high and has less time delay. If there would be height information available, then it 
would be possible to estimate the propagation speed of oscillations from one height to another. The radiance 
va riation of several line s with respect to solar-Y has been plotted in Figure|4]at solar-X~ -72". As described 



m 



Q'Sheaet all (|2006r) . the difference in the peak positions provides an estimate of the differences in the 
formation heights between different lines at that particular time and condition. Hence, using the correlation 
technique, the time delay between two lines can be obtained and using the limb brightening technique, the 
formation height difference can be estimated for a line pair. The results obtained have been summarized 
in Table [2l Due to the relatively poor temporal resolution, uncertainties in time delay measurements are 
large. Hence, for the on-disk bright region, results obtained here can only be used to infer that waves are 
propagating from lower to higher heights in the solar atmosphere. 



Table 2: Linear correlation coefficients between oscillations in different line-pairs corresponding to the on- 
disk bright region. 



Line- 


log Tmax 


Correlation 


Time delay ^ 


Height diff.2 


pair 


(K) 


co-efficient 


(s) 


(km) 


He ii/Fe xii 


4.9/6.1 


0.242954 


188 


4509 


He u/Fe xi 


4.9/6.1 


0.250901 


94 


- 


He ii/Fe x 


4.9/6.0 


0.322859 


47 


5628 


He u/Fe vm 


4.9/5.6 


0.307703 


47 


3180 


He ii/Si VII 


4.9/5.8 


0.564139 


47 


5363 


He ii/Mg VII 


4.9/5.8 


0.640592 


47 


4104 


He ii/Mg VI 


4.9/5.6 


0.496354 


47 


3654 


S v/Ne VIII 


5.2/5.8 


0.504659 


36 


3647 


iv/Ne VIII 


5.2/5.8 


0.477825 


54 


2917 



^Limited by time resolution defind by effective cadence of respective instruments. ^Limited by spatial 
resolution of respective instruments, w 715 km. 



A similar correlation analysis is applied to the offlimb inter-plume region at Xw -72" using data from 
the same line, Fe xii, but at different heights with respect to solar- Y~ 1000". The results aie plotted in the 
right panel of Figure [13] The time resolution is about « 47 s, the EIS cadence. Also in this case, the level 
of correlations and time delays are in inverse proportion as expected. The measured time delays are plotted 
against solar- Y in left panel of Figure [14] The continuous line corresponds to a second order polynomial fit 
applied to the data points (as marked by asterisks). The error bar on these time delays is obtained from the 
Half Width Half Maximum (HWHM) of the particular correlation plot. The dotted line corresponds to the 
fit to the slanted radiance ridges in the x-t maps (white lines of Figure [6]). The figure indicates that the travel 
time is decreasing with height indicating an acceleration, as was seen from the x-t map of Fe xii in Figure [6] 
This Figure provides an independent estimate of the acceleration. 

The Alfven wave speed through the qu asi-static corona can be calculated from the expression Va = 



B/ ^JAnp. We use the density profile given by Teriaca et al. j 2003 ) and take into account the superradial fall 



of the magnetic field with height according to JKopp & Holzen (Il976r) . with a base magnetic field of « 0.65 G. 



For comparison purpose, the Alfven wave speed is also calculated by assuming a constant magnetic field of 
w 0.65 G with height. The measured time delays at different Solar-Y in inter-plume region are then compared 
with the travel time for a theoretical Alfven mode. Th e Alfvenic time delays obtained assuming a magnetic 
field constant with height and expanding according to iKo pp & Holzen (11976 ). are plotted, respectively, as a 



dashed and a dot-dashed line in the left panel of Figure [14] To estimate the propagation speed, we plot right 
panel of Figure [14] which shows the variation of propagation speed with respect to the height in the solar 
atmosphere. The speed is calculated by the time derivative of the fit to the measured time delays shown in 
the left panel of Figure [14] and is plotted as a continuous line. This speed is compared with the theoretically 
calculated propagation speed of Alfven modes, plotted as dot-dashed and dashed line, respectively, for the 
expanding and constant field case. It can be seen that the measured propagation speed is roughly consistent 
with being Alfvenic if we assume a field of 0.65 G at the base. From the figure it can be seen that near the 
limb and off-limb, the speed of propagation is about 130 km s~', increasing to more than 220 km s~' far 
off-limb, close to the speeds obtained from the x-t map in Figure [6] It also appears that above a certain height 
(Y ~ 1080") the speed increases more rapidly. These results might indicate that physics of the propagation 
might also change at these heights. 



4. Discussion 

Our results show the propagation of disturbances from an on-disk region to the off-limb corona within 
inter-plume region. These disturbances appear to originate from an on-disk bright location, presumably 
the footpoint of a coronal funnel, around solar-Y~ 967" (see Figure [5]). No signature of propagation is 
visible below solar-Yw 967". The propagation speed, as measured from the x-t map (distance-time map) of 
Ne vm and Fe xii line radiance, increases with height. Outward propagating disturbances are also recorded 
in a plume structure within the EIS field of view, but the acceleration is almost zero. Furthermore, the 
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Fig. 13. — Left: Correlation coefficients versus time delay between the line-pairs at the on-disk bright 
location (X ~ -72", Y« 961"). The maximum correlation coefficient for a fixed line-pair provides a 
measure of the travel time. Right: Correlation coefficients versus time delay in inter-plume region for time 
series in the EIS Fe xii line. The correlation coefficients are calculated at different heights with respect to 
Y« 1000". 



disturbances become diffuse far off-limb as seen from the EIS Fe xu x-t map. From wavelet analysis these 
disturbances have a periodicity in the range of 15 min to 20 min in both regions and the periodicity is 
seen almost over the whole duration of observation, although the power is increasing and decreasing with 
time. A correlation analysis of the light curves of different spectral line-pairs in the on-disk bright region 
reveals larger time-delays for line pairs with the larger difference in formation height (inferred from the limb 
brightening curves), indicating upward propagation. This upward propagation indicates that these waves are 
generated somewhere lower in the atmosphere, probably in the chromosphere, since these are also seen 
in He ii and then propagate upward towards the off-limb region. In the off-limb region, cross-correlation 
between light curves from the same spectral Une (Fe xu) but obtained at different heights above the limb 
shows time-delays indicating outward propagation, in agreement with the results from the analysis of the x-t 
maps. 

The propagating disturbances, clearly visible in the inter-plume region in the Ne viii radiance x-t map, 
are not visible in the Ne vni Doppler shift and width x-t maps. However, the wavelet analysis shows that 
both Doppler shift and width oscillations are present in the on-disk bright region and near off-limb with 
approximately the same periodicity as seen in radiance (see Figures [H |9l &[TOl). and they are also mostly in 
phase. The presence of oscillations in radiance and both resolved (Doppler shift) and unresolved (Doppler 
width) velocities with approximately the same period, is evidence of propagating waves with at least a 
compressional component in the inter-plume lanes. The measured speed of propagation of these waves from 
the Ne vm line is (25 ± 1.3) km s"', (38 ± 4.5) km s~' and (130 + 51) km s"' in the on-disk region, 
near the limb and off-limb respectively. Fe xu data shows that the propagation speed further increases to 
(330 ± 140) km s ' in the far off-limb region of inter-plume. In the plume region, instead, the observed 
off-limb propagation speed increases from (135 ± 18) km s~' to only (165 ± 43) km s~' far off-limb. 
Beyond this point the radiance disturbances become diffuse. The increase in propagation speed is small and 
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Fig. 14. — Left: Variation of travel time with height in inter-plume region. The asterisks represent the 
measured time delays from right panel of Figure [13] while the continuous line corresponds to a second order 
polynomial fit applied to the data points. The error bars on these time delays are obtained from the Half 
Width Half Maximum (HWHM) of the respective correlation peaks. The dotted line corresponds to the fit 
to the slanted radiance ridges in the x-t maps (white lines of Figure O. The dashed and the dot-dashed 
lines are the theoretically predicted Alfvenic time delays o btained assuming a magnetic field constant with 
height and expanding according to Kop p & Holzen (119761) . respectively. The change in the slope indicates 
acceleration. Right: Propagation speed with height in inter-plume region. The continuous line is obtained 
after taking the time derivative of the fit to the data points in the left panel. The dashed and the dot-dashed 
lines are the theoretically p redicted Alfven speeds obtained assuming a magnetic field constant with height 
and expanding according to lKopp & Holzen (119761) . respectively. 



21 



the acceleration is negligible within the given uncertainties. 

The propagation speed becomes supersonic (> C., ~ 170 km s"\ for Fe xii line formation temperature) 
far off-limb in the inter-plume region. Moreover, near the limb region, the Ne vni and Fe xu lines, which are 
in phase, show nearly equal propagation speeds for some overlapping region despite having quite different 
formation temperatures, suggesting that these propagating disturbances are temperature independent. This, 
together with the presence of oscillations in Doppler width and shift, may suggest that these waves are 
Alfvenic in nature. Pure Alfven waves do not cause any density perturbation and thus would not cause 
any radiance fluctuation. However, in the real case of waves propagating in a density stratified atmosphere, 
non-linear effect may also cause small densi ty fluctuations leading to radiance changes of a few percent 
(ITu & MarschI 1 19951 : iKaghashvih et all 120091). Moreover, oscillations in the obser ved line widths can be 



caused by torsional Alfven waves (|Zaqarashvilil 12003 



Van Doorsselaere et al. 



20081). Hence, interpretation 



of these propagating disturbances in terms of Alfven waves appears quite reasonable. 

Furthermore, observed radiance oscill ations can also be e xplained as a line-of-sight effect of entirely 
incompressible MHD waves as described inl Cooper et al.l(l2003b . In this model, when observed at an angle 6 
to the direction of propagation, the wave-induced deformation in a coronal loop causes brightness variations. 
This is because the amount of optically thin emitting plasma along the line-of-sight changes as a function of 
time. Thus, the radiance variations can be produced even by entirely incompressible MHD waves. 

Conversely, the measured propagation speeds are also consistent with the fast magnetoacoustic mode 
of propagation within the error bars of the propagation speeds in Ne vni and Fe xii lines and explains the 
observed radiance oscillations due to its compressible nature. Hence, interpretation of these propagating 
disturbances in terms of fast magnetoacoustic waves also appears reasonable. 

It should be noted that it is always the apparent propagation speed in the plane of the sky that is 
measured. The structures carrying the waves most likely form an angle with the plane of the sky and, thus, 
these propagation speeds are always a lower limit. However, in the case of the inter-plume region, we trace 
the origin of the wave to an on-disk bright region located at Xw -72" and Yw 967". Assuming roughly 
radial propagation, this would imply an angle of « 9° with the plane of sky, leading to a difference between 
apparent and real propagation speed of less than 2 %, a negligible quantity. In the case of the plume, a 
careful inspection of Figure |2] shows that the plume is rooted very close to the visible limb and is lying close 
to the plane of the sky. This suggests that in this case the difference between apparent and real propagation 
speed is likely very small. 

The k ind of oscillations which are presented here are very similar to other already reported in the liter- 
ature (e.g. 



DeForest & Gurman^ll998 



Q'Sheaetal. 



Ofman et al. 



20001: Baneriee et al.ll2000al . l200ll : Morgan et al.ll2004l : 



2006.. .2007. : .Gu pta et al.ll2009. : .Baneriee et al...2009b.) in on-disk and off-limb regions of coro 



nal holes observed with different instruments. iBanerjee et all (1200 ih have reported osc i llation s in on-disk 
network regions and in inter-plume regions similar to those reported here. lOupta et all (120091) have calcu- 
lated the speed of propagation in the on-disk bright region using the statistical technique and the measured 
speed were of the sam e order as reported her e. All these observations indicate the presence of propagating 
MHD waves. Recently, iMcIntosh et all (|2010h detected propagating features in polar plumes using STEREO 
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observations. These authors have interpreted these features in terms of high speed jets of plasma traveling 
along the structures which repeat quasi-periodically, with repeat-times ranging from 5 min to 25 min. This 
would be contrary to the widely held interpretation that this observational phenomenon is due to compressive 
waves. One should also consider Doppler dimming results from UVCS and SUMER aboard SoHO, which 
shows that velocities above 100 km s~^ are re ached only above 1.5 R/R© in either plumes or inter-plumes 



(iTeriaca et al.ll2003l: iGabriel et al.ll2003L l2005h . 



To our knowledge, this is the first time that a signature of accelerating Alfvenic waves or fast magne- 
toacoustic waves originating in an on-disk bright region has been observed in near and far off-limb regions 
within 1.2 R/R©. In inter-plume region, the disturbance is seen propagating up along the whole EIS slot 
length whereas in the plume region, it becomes diffuse far off-limb. In inter-plume region the wave prop- 
agates farther in the corona with high acceleration whereas in plume region it may have been dissipated in 
the off-limb region. However, the lack of signature at greater heights in the plume may be simply due to 
merging with the background signal. This suggests that inter-plume regions may be the prefer red channel 
for the acceleration of the fast solar wind. This conclusion is in agreement with earlier reports (|Wang et al. 
I997I : IPatsourakos & Viajboool : loiordano et aPboOol : baneriee et aPboOll : iTeriaca et allboOBh . 



5. Conclusion 



Analysis of Ne vra and Fe xii radiance x-t maps reveal the presence of outward propagating radiance 
disturbances in the off-limb and neai" off-limb region of inter-plume with periodicities of about 15 min to 
20 min. From the SUMER Ne viii line radiance x-t map, one can infer that the waves originate from a 
bright location (presumably the footpoint of a coronal funnel) and propagate towards the limb with a speed 
(25 ± 1.3) km s"^ Aroundthelimb the speed has increased to (38 ± 4.5) km s"^ reaching (130 ± 51) km s~^ 
off-limb. Further far off-limb, the speed of the propagation becomes (330 ± 140) km s'^ as seen in the 
EIS Fe xn line. Similar propagating disturbances are also seen in the plume region but with negligible 
acceleration, if any. The waves are not visible far off-limb, suggesting that they may be dissipated or, more 
simply, merge into the background. The waves as recorded in the inter-plume regions are either Alfvenic 
or fa st magnetoacou stic in nature whereas the one seen in plumes are more likely slow magnetoacoustic 
type. ITu et all (120051 ') have conjectured that the solar wind outflow is launched by reconnection at network 
boundaries between open flux lines and infra-network closed loops. The infra-network closed loops are 
pushed by supergranular convection towards the network triggering reconnection. This scenario is consistent 
with our identification of the origin of the propagating disturbances in inter-plume region with an on-disk 
bright region. These results support the view that the inter-plume regions are the preferred channel for the 
acceleration of the fast solar wind. 
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